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Abstract 



We analyze the effect of adding quenched disorder along a defect line in 

the 2D conformal minimal models using replicas. The disorder is realized 

rj-^ ' by a random applied magnetic field in the Ising model, by fluctuations 

f-H , in the ferromagnetic bond coupling in the Tricritical Ising model and 

Q ' Tricritical Three-state Potts model (the (j>i2 operator), etc.. We find that 

O ' for the Ising model, the defect renormalizes to two decoupled half-planes 

without disorder, but that for all other models, the defect renormalizes to 
a disorder-dominated fixed point. Its critical properties are studied with 
an expansion in e oc 1/m for the m*'^ Virasoro minimal model. The decay 
exponents Xat = f (I - j^^^ + 0{-^f) of the iV* moment of the 
OO , two-point function of <j!>i2 along the defect are obtained to 2-loop order, 

exhibiting multifractal behavior. This leads to a typical decay exponent 
^typ = ^(1 + -, — ^ryy + 0{ ^,-^ )'^). Oue-poiut fuuctious are seen to have a 
non-self-averaging amplitude. The boundary entropy is larger than that 
of the pure system by order 1/m'^. 

As a byproduct of our calculations, we also obtain to 2-loop order the 
j^ . exponent Xiv =7V(I- g|^(3iV-4)(g-2)2-FO(g-2)^) of the Af*'' moment 

of the energy operator in the q-state Potts model with bulk bond disorder. 



1 Introduction 

Conformal symmetry tends to emerge in pure (homogeneous and rotationally 
invariant) 2-D Statistical Mechanics models at their critical points. This high de- 
gree of symmetry severely constrains these theories, so that these critical points 
are well understood. Many models, including the Ising model, Tricritical Ising 
model, and Tricritical Three-state Potts model, are part of a class of conforrnal 
field theories known as Virasoro minimal models H, ||, y| . The Ising model 
Tricritical Ising model [^ y, ^, and Tricritical Three-state Potts model |£ 
have all been realized experimentally in adsorbed monolayer systems. 

Because physical systems always have impurities, it is important to consider 
the effect of quenched disorder on the critical behavior of these theories. When 
disorder is added to the bulk, random fields are usually relevant, but random 



bonds may or may not be (see e.g. Ref.'s y, 10 ). It is possible to show, 
for certain systems, that with the addition of disorder the system renormalizes 
into an infra-red fixed point. In fact, a rigorous theorem shows that when 
disorder in the order parameter is added to a system undergoing a first-order 
phase transition, the latent heat vanishes O], and a 2"^^ order transition can 
be expected. 

However, one generally does not know a priori whether the new critical point 
is disorder-dominated. A number of studies have reported cases in which the 
addition of quenched disorder to a pure 2-D model at its critical point resulted 
in another pure (that is, non-random) critical model pa O, nj, na, O, nTj. 
Other studies found quenched disorder to result in new disorder-dominated fixed 
points Pq, 19l 00, Ell 03]. One can often see that a critical point is disorder- 
dominated by showing that various universal quantities are not self-averaging 
(that measurements for a specific 'typical' sample may differ substantially from 
those on an average one). One particularly interesting manifestation of this is 
multifractal behavior, which occurs when an infinite hierarchy of independent 
scaling dimensions are associated with a single operator ||2^, gj, g5|. (See for 
example [^ for a relevant discussion.) 

All the above studies have focused on the effects of adding quenched disorder 
to the bulk of a system. However, it is also interesting to consider the case where 
the 2-D model has a defect along which impurities have clustered. In this paper 
we consider the effects of adding quenched disorder only along a defect line, in 
each of the minimal models; these models are labelled by an index to, ?n > 3. 
Each m represents a different model : the Ising model {m = 3), Tricritical Ising 
model {m — 4), Tetracritical Ising model {m — 5), Tricritical Three-state Potts 
model (to = 6), etc. . . In Section we introduce our defect model, adding 
quenched disorder in the coupling to 0i2 (an operator in the Kac Table [|l[) 
using replicas. The physical meaning of this disorder varies from model to model, 
representing a random magnetic field for the Ising model, but fluctuations in the 
ferromagnetic coupling (or chemical potential) in the Tricritical Ising model and 
Tricritical Three-state Potts model. We find that the Ising model renormalizes 
to a non-random model (consisting of decoupled half-planes with free boundary 
conditions and no random magnetic field), while all other models renormalize 



to disorder-dominated fixed points. Just as with bulk disorder, disorder on the 
defect can result in either a pure or disorder-dominated fixed point. 

In Section we calculate the renormalization group equation of the strength 
of the disorder A, to 2-loop order by minimal subtraction p^, EtI. We find 
that arbitrarily weak disorder grows, and flows to a new fixed point at which 
A is of order 1/m. This justifies a 1/m expansion, where critical quantities are 
calculated perturbatively in 1/m. The "boundary entropy" [Esj at the random 
fixed point is found to be 0{\/m?), and larger than the entropy of the pure 
fixed point. 

In Section ^ the same scheme is used to find the moments of correlation 
functions of the operator 0i2 along the defect. (Technical details associated 
with the irreducible representations of the symmetric group |^ are delegated 
to Appendices ^ and y.) It is found that the moments fall off as a sum of power 
laws, the dominant term decaying as 



< 0i2(xi)0i2(a;2) >'^ (x| XI - X2 r'^" (1.1) 



with 



/V / 1 '^ S \ 

^^ = ^ 1 - -AiN 4)(^)^ + 0{^f (1.2) 

2\ 4 m + 1 ni + 1 J 

Xm < NXi, so the operator <j>i2 exhibits multifractal behavior. This leads to a 
typical decay exponent [E9| : 



Xt.pic.4(l + (^)^ + C^(^)^) (1.3) 

In Section we calculate one-point functions of the same operator off the 
defect line. The universal (normalized) amplitudes of the moments are found 
to be non-self-averaging, wheras the power law is self-averaging. 

In Section o, the Ising model, which requires special considerations due to 
the presence of an additional marginal (boundary) operator, is analyzed. It 
is argued that when a random magnetic field is added along the defect line, 
the system renormalizes to two decoupled half-plane Ising models with free 
boundary conditions and no disorder. 

Wc finally note that the manipulations needed to get the exponents in 
Eq.(n]3) arc similar to those needed to get the decay exponents of moments 
of two-point functions of e, the energy operator, for the g-state Potts model 
with bulk disorder in the bond strength. This model has been analyzed else- 
where by expanding about q — 2 (the Ising Model) (2l|, |2^ , P9|, P7|, and as a 



byproduct of our calculations here, we find in Section ^ Eq.( 4.11 - 4.1S ), that 



<s{xi)s{x2)>^ 'x\xi-X2\-'''" (1.4) 



where 



Xn^N{1- ^{3N - A){q ~ 2f + 0{q - 2f ) (1.5) 

Results of numerical transfer matrix calculations for the random bond q- 



state Potts model can be found in |22 



2 The Defect Model 

We start with a Virasoro minimal conformal field theory |^, g, ^ labelled by an 
integer ttt,, ttt, > 3, and perturb it along a defect line. The perturbed action is 

/CX) 
dx/i(x)0i2(2;,y = O) (2.1) 

-oo 

Sm, is the action of the unperturbed conformal field theory. h(x) is a random 
coupling, is picked from a Gaussian probability distribution with zero mean and 
variance Aq, and is uncorrelated along the defect: 



/i(a;)=0, h[x)h{x')^2l\Q8{x-x') (2.2) 

The overbar indicates the disorder average. The operator 0i2 is located at 
position (p, g) = (1, 2) in the Kac Table |]l|, and exists for any minimal model. 
For the Ising model (ra — 3) it is the spin operator, while for the Tricritical 
Ising model {m = 4) and Tricritical Three-State Potts model {m — 6) it is the 
energy operator. The scaling dimensions 2hpn (twice the conformal weight) of 
operators (t)pq in minimal models are known 1I| : 



[{m + l)p — mq\'^ — 1 
2m(m + 1) 



^K. = " .:^L :i > (2-3) 



which gives 



2hi2 = --^-^ r (2.4) 

2 2(m+l) ^ ^ 



The replicated |Q effective action is 

grepUca ^ ^ ^a ^ ^ / ^^ ^ ^a^^^^ {x , y ^ 0) , (2.5) 

where n is the number of replicas, and we take n — > at the end of the calcu- 
lation. 



We have ignored higher cumulants of the probabihty distribution of h{x), 
because power counting shows that they are irrelevantly in the R.G. sense, for 
large to, and we will be expanding about large to. 

We have also dropped the terms with a = /3 in Eq.( |2.5[ ). The terms with 
a = P produce a non-random (/)i3 by the conformal fusion rules [Q. Upon 
renormalizing they will generate other non-random (j^iq with g > 3. However, 
from Eq.( p.4[ ) we can see that all these terms ar e irr elevant except when q = 
TO = 3. So when to ^ 3 the perturbation in Eq.(2.5) is the only relevant one. 



The following analysis in Sections |^ ^ will assume this, and will thus only hold 
for m y^ 3. This will not be too restrictive, since all quantitites will be calculated 
by expansion in e = — xr, and will thus be based on large to,. But if to = 3 
(the Ising model), then ^13 (the energy operator) is marginal, and we need to 
include the effects of a constant coupling to 0i3. This is done in Section ^. 

3 Renormalization of the disorder strength A 

We calculate the renormalization group equation for A to 2-loop order by min- 
imal subtraction pl|, p^. To calculate the renomalization of Aq, we want to 
know, given a microscopic disorder strength Aq, what effective disorder strength 
A(r) is seen on large length scales r. In a region of size r, we expand out 
expjS'^P'*'^''} in powers of Aq. (Aq)'', with p > 0, will come with p disorder 
operators at various points in the region of size r. p = 2, 3 or more disorder 
operators may look, using repeated operator product expansions, like a single 
disorder operator on larger length scales (or size r), and thus create a new 
effective disorder strength A(r). This can be represented schematically as : 

/ \ 2 

dx J2 <t>l2<t>U^) + o ^0 / dxY, '^?2'^?2(^) + 

+ - Ao / dxY,<t>'iX^{x)\ +...^A dxJ2<l>i2^i2i^) (3.1) 

This analysis will of course generate numerous terms besides the disorder 
operator. However, as noted above, these terms are all irrelevant in the RG 
sense, so we will not calculate these terms. 

For each power of Aq, the integrals generated above are regulated at short 
distances by analytic continuation in e = ^.^ , and at large distances by an 

^ For the lower values of m = 4 or 5, the 4"^ order cumulants are respectively relevant 
and marginal, but the 2"^^ order cumulant is more relevant, and presumably will be physically 
dominant. While the 4*'^ order cumulant is not irrelevant at the unperturbed fixed point for 
m = 4 or m = 5, it is irrelevant for large m, and we expect it to be irrelevant for m = 4 or 
5 at the new disordered fixed point. This would be similar to the Wilson-Fisher fixed point, 
where the 4>^ operator becomes relevant at the Gaussian fixed point below three dimensions, 
wheras it is in fact irrelevant at the new fixed points obtained by an expansion about four 
dimensions. 



infrared cutoff r. The calculation, which uses the method of [ETI, is done in 
Appendix H, where we obtain 



A(r) =r'Ao+4(n-2) — Ag-4(n-2) 



2-4(n-2)- 
e 



-Al + OiAt) (3.2) 



When we calculate the beta function by taking a derivative with respect to 
log(r), we find that the poles in e cancel, as they must for any physical quantity. 
The result is 

/3(A) = — ^^ = eA - 8A2 + 32A3 + 0(A4) (3.3) 

d(log(r)) 

where we have taken the replica limit n — > (and have ceased writing the r 
dependence of A(r)). The RG flows take the unperturbed theory to a new 
infrared fixed point with disorder A*. Solving for A» by putting /3(A*) = 
gives 

A,44+0(e^), ie^^^) (3.4) 

The new fixed point is at a distance of order e from the unperturbed theory, and 
so we can calculate physical quantities by expanding in powers of e, which is 
small for minimal models with large index m. This is analogous to the Wilson- 
Fisher epsilon expansion in d = 4— e dimensions. We will see below that the new 
fixed point is disorder-dominated. Because the analysis is based on expansion in 
powers of e = ^.-^ , we really only show that there is a new disorder-dominated 
fixed point for large m - but we expect no qualitative changes for lower m, such 
as m = 4. However, we again note that to = 3, the Ising model, is qualitatively 
different, and will be treated separately. 

The 'boundary entropy' (the universal constant independent of the system 
size appearing in the disorder-averaged free energy) associated with the defect 
line, calculated as in Ref.[g^, is found to be 

Note that the entropy has increased from that of the unperturbed system (where 
it vanishes). This is to be contrasted with the case of a pure (non-random 
and unitary) system, where the entropy ||28| is expected to only decrease upon 
renormalization. 



4 Scaling Dimensions of Moments of 



h2 



We now look at < 012(2:1)012(2:2) >", the disorder-averaged N^^ moment of the 
2-point function for points xi and X2 which lie both near the defect and far 
from each other (Figure 1). We want to see how these moments fall off at large 
distances. As formulated in replicas we have 



N N 



< 012(xi)012(a:2) >'^ = ( n <^?2(^l) n <^12(^2) ) (4.1) 



/3=1 



tN 



As explained in p9| |, we cannot simply calculate the dimension of 110= i '/'?2; 
because this operator is not multiplicatively renormalizable. Instead, it is a sum 
of independent scaling operators with different scaling dimensions. In the limit 
where the pure and disordered fixed points collide {m = oo), we have numerous 
operators with the same dimension. For example, looking at TV = 2, 0j2<^i2 
and (t>i24'i2 have the same scaling dimension and are equally good operators 
(note that the replica limit n — » is not taken until after all calculations are 
completed). In general, for the A''*'' moment, and with n replicas, we have (^) 
operators with the same dimension at m = cxd ; because they all have the same 
scaling dimension at ?7i = cx3, any linear combination of them also has this 
scaling dimension at ?n = cx). 

However, when we move to ?7i < oo, only appropriate linear combinations of 
these operators will have well-defined scaling dimensions. The entity Jla=i "^12 



which appears in Eq.(l.l) will contain all of these scaling operators, so that the 
average < *i2(»:i)0i2(=:2) >o" will decay as a mixture of power laws, and at large 
distances will be dominated by the scaling operator with the lowest dimension. 
We thus need to calculate the dimension of each of these scaling operators. The 
appropriate multiplicatively renormalizable operators transform in irreducible 
representations of the symmetric group p9| : 

l<Qi<(ri-J\/) 

where < M < TV. 

To calculate the scaling dimensions of On Mm we will add the term 
^NMnfl ^_ dxONMn{x) to the action. As in section y, we calculate the reno- 
malization of Atva/a^.o to IS.NMni'i') — Z n Mnir) ^n Mn.o on length scales of size 
r, by expanding the action in powers of Aq: 



^NMnfi ( / dxONMn{x) 1 + 

+ AArM«,oAo f / dxONMn{x)\ / dx' ^ 0"2 012 (2^' 

\J-oo J \J-oo ^_,p 

+ -ANMn,oAl( I dxONMnix)] f dx' J^ ^12(l^l2ix') 




ZNMn{r)/^NMn,0 ( / dxONMn{x) 1 (4.3) 



In Appendix H we check that to 2-loop order we don't need to worry about 
mixing with other operators. ZNMnir) is calculated with the same type of in- 
tegrals as in the last section (they are again regulated at short distances by 
analytic continuation in e and at large distances by an infrared cutoff r) . How- 
ever, technical combinatorial complexities arise in counting the number of con- 
tractions associated with various irreducible representations of the symmetric 
group. They are delegated to Appendices and 0. The result is 



j.2e 



ZNMnir) = 1 + 2bNMn-^0-'i[N{n-N) + {N-l)bNMr. 

+ {2{bNMn)^ + 4(n - 2)bNMn) ^ Ag + ©(Ag) (4.4) 

where we have defined 

bNMn = 2 ((iV - M)n -N^ + M{M - 1)) (4.5) 

iiNMn is a solely combinatorial factor which arises from counting the number 



of replica contractions associated with Onmu- We then use Eq.(3^) to rewrite 
the series in Aq as a series in A, yielding 7ArMn(A). 

(.. _ d(log(ZjVM„(r))) 
d(log(r)) 

= 26wM«A - 8(iV(n ~ N) + [N - l)bNMn)^'' + 0(A3)(4.6) 

Note that the poles in e again cancel. To get the scaling dimensions of Onmu 
at the disordered fixed point, we take n — > and A — » A», getting 

7jVAf (A,) = ^bNMo + ^{N^ -{N- 2)bNMo) + 0{e^) (4.7) 



In the unperturbed theory, for all M, < Onm7^(,xi)OnmA^2) >" will fall off as 
I a;i - a;2 j-^^'^i^. But with the defect it will fall off as | xi - X2 |-2^"Af , where 
2Xmm — 4:Nhi2 ~ 27ArM(A»). Because na=i <^12 i^ ^ linear combination of the 
scaling operators Onmu, the moment < 4,i2(=i=i)4>i2(^2) >™ will be a sum of terms 
decaying with powers X^m- For | a;i — X2 | large, this will be dominated by the 
smallest power, and it is easy to see that this is X]\j = XjsfN, giving our main 
result : 



< 0i2(xi)0i2(a;2) >" oc| XI - X2 \-'^" (4.8) 

N f2 

XN^-{l~-{3N-4) + 0{e')) (4.9) 

Xn < NXi, so the system is not self-averaging. Instead, we have an infinite 
number of independent scaling dimensions all associated with the single operator 



(f)i2. Note that X^^/Ni > Xn^/N2 for N2 > Ni, as required by convexity. The 
result above also yields the typical |29] exponent: 



X, 



typical 



^<' 



0(6^)) 



(4.10) 



The combinatorial problems encountered here are the same as those for the 
q-state Potts model with disorder in the bulk ferromagnetic couplings. The 
specific integrals are different but have already been done in ||27| . The only 
new difficulty is that to two loop order, the operators Onmu are not longer 
always multiplicatively renormalizable, but instead mix with other descendent 
operators. However, luckily, the leading and subleading operators, Onnu and 
ON,N-i,m remain multiplicatively renormalizable (see Appendix]^ for details). 
So the combinatorics in Appendices ^ and ^ also give the 2-loop result for 
the decay exponent of the A^'^ moment of the two-point function of the energy 
operator e, in a Potts model with bulk bond disorder. We have: 



< £{xi)e{x2) >^ oc| xi - a;2 I 



-2Xa 



Xn = N{1- —{3N -4){q- 2)2 + 0{q - 2)^ 



This also gives us the typical decay exponent : 



(4.11) 
(4.12) 



^typical = 1 + -^(g - 2)2 + 0(g - 2)3 



(4.13) 



5 Moments of the One-point Function of 



^12 



We now calcuate the disorder-averaged moments of the one-point function of 
(j)i2, evaluated at a distance y from the defect Hue (Figure 2). The method is 
similar to that in |3^| . We no longer need to worry about the various irreducible 
representations of the symmetric group, because the one point function oiONMn 
vanishes by symmetry except when M = 0. The disorder-averaged odd moments 
of 012 vanish because the disorder-averaged system is symmetric under 0i2 —^ 
—(f>i2- The even moments are calculated perturbatively : 



612(2; = 0,2;)) 



2N 



2N 



Yl^i2io,y) 



2N 



n'^?2(0,y) 



„'^«/_°°J^Eo^0'^?2'?if2 



,-(1- 



^-^^"E^(Aoj 



^N 



■27V 



.Q = l 



0, cutoff a 



n'^?2(o,i) 



n / "^^j 



OO -. 



n ^fi-.--.i-^ 



(fe<0=i 



i=M 



y 



' 

(5.1) 



In the 2"*^ line we have introduced a cutoff a to regulate the short-range 
divergences. In the 3''^ line we have expanded out the effect of the defect pertur- 
batively in Aq, and used the conformal symmetry to rescale each 
expectation value by y. In the lowest order term, Ij^ , to get a nonzero expec- 
tation value, the 2N operators on the defect must lie in different replicas - so 
nothing special happens when two defect terms get close together, and we can 
drop the cutoff a. 



/r(-) = (2iV)! 



y 



dx 



>(2^^ + l) 



l-e 



N 



= (2iV)!7r^(l + 0(-,e)) 
V 



(5.2) 



We can now use the operator product expansion (OPE) 



Y,^i2ix + 5,y^mi2i^ + 5.y = ^) 



a#/3 



^ 0?;(x,y = 0)0^2(3^:2/ = 0) 



a'^13' 



^^"^ ^^E'^?2(a:,2/ = O)0f2(x,y = O) for 5^0 (5.3) 



(51- 



a^l3 



to get the leading poles in e for the higher-order terms, by taking derivatives 
of l\ (-) with respect to -. When we take derivatives of the step function 
9{\ Xi — Xj I — -), we bring two operators close together, and so can use the 
OPE. We get 



^i-^=-4^(^-l)(n-2)(J)-l+^/f)(J) for ^>N (5.4) 



AN) 



(N) 



We already have 7]^ , so we can now repeatedly integrate to get l!> ' ior i > N. 
/j- must be finite as e ^ for a ^ 0, because with a finite cutoff a, there are no 
ultraviolet singularities in the integrals for I^ - this requirement determines 
all the constants of integration. So, for example, the leading term in -?^jv+i(~) i^ 



AN) 



AN) 



I'nM-) - -4iV(iV+l)(n-2)/]^^^(0) 



1 ,a 

-[ — ) + constant 

e y 



(l + 0(-)) 

y 



-ANiN + l)(n - 2)/W (0) [ ^^^ ) (1 + 0{^)) 



10 



e^O 



^4Ar(iV + l)(n-2)/(f)(0)log(-) 



(5.5) 



(Note that this does diverge as a — > 0.) More generally, by repeatedly integrating 
Eq.(^.4|) and using the requirement that all terms be finite as e ^ for a ^ 0, 
wc find 



'.'"'<?) -ir^iiM^(-^<"-) --9 



i-N 



4^0), ^>N (5.6) 



Now taking the leading divergences of each term in Eq.(5.1) gives 



(012(2/)) 



2N 



m(N- 1)1 ^{i-N)r ^ ' ^y"^ °' 



^N (0) ^(l-e)JV 

iV! ^ 



Aoy^ 



N 



l + 4(n-2)log(|)y'Ao^ 



(5.7) 



In the front we have extracted a constant and the expected power-law de- 
pendence y~^^~'^>^ = y-'^hi2N ^ The remaining terms are written as a function 
F[a, y, Aq] of the large distance y and the microscopic variables a and Aq. How- 
ever, we know from the Callen-Symanzik equation that this can be rewritten in 
terms of a single scaling function dependent only on the renomalized coupling 
for length scales of order y, A(log(|)) Q. Exphcitly, 



F[a, y, Ao] = F[ae', y, A(^)] = F[y, y, A(log(^))] = G[A(log(^))] (5.8) 



We rewrite Aq in terms of A(log(|)) by integrating Eq.( |3.3|) to 1*** order, getting 



Aoy^ = A(log(^)) -f 81og(^)(A(log(^)))2 + 0(A(log(^)))3 + . . . (5.9) 
a a a a 



Substituting this into Eq.(5.7), we indeed get that the leading divergences at all 
orders of perturbation theory sum up to give a function which depends only on 
A(log(|)) : 



2iv_ {'2N)\ ^^A(log(f)) 



^{y)Y" = 



i/2/ii 



N 



(</>i2(2;))''^+' = 



(5.10) 



We see that the amplitudes of one-point functions are not self-averaging : 



(0i2(y)) 7^ [{4>i2{y)) ]^- That is, the average of the N^^ power is different 
than the N^^^ power of the average. The amplitude ratios are universal properties 
of the random defect fixed point. 

Also note that ^^ log [ J^mHW^' 1 < i^ log [ JMyTf^' 1 for ^i < A^a, 
as required by convexity. 



JVa 
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6 Ising Model 

We now look at the defect introduced in Eq.( |2.1| ) for the special case of the 
Ising model. As noted at the end of Section 0, the analysis in the three sections 



above fails when m — 3, because along with the perturbation in Eq.(2^), we 
generate the marginal non-random operator 0i3. The minimal model with m = 
3 corresponds to the Ising model, and in this model 0i2 is the spin operator 
(a) and <j)i3 is the energy operator (e) . So in this case our defect consists of a 
random applied magnetic field and a non-random energy operator. Our action 
is 



/"OC 

D = Olsing + / "2; 



A Y^ cr„(x, O)fT0(x, 0) + A ^ e„(x, 0) 

a^P a=l 



(6.1) 



In the lattice formulation of the Ising model, perturbing with the energy 
operator is the same as changing the bond strength. Defects where the bond 
strength is changed along a single line have been studied and solved exactly 
by Bariev and McCoy et. al. [^ Q. They dealt with cases where the bond 
strength was changed only in the bonds perpendicular to the defect (the ladder 
geometry - Figure 3), or only in the bonds parallel to the defect (the chain 
geometry - Figure 4). Our defect, for the Ising model, is thus a perturbation 
with a random magnetic field, of these exactly solved defects. The scaling 
dimension of the spin operator along the defect is Xa- — |(?(A)^, where 



2 / /2 — 1 \ 

5(A) = — tan^^ — -— — for the ladder geometry, (6.2) 

TT \ tanh(Ac + A) / 



and 



g[\) = —cos ^ (tanh(2A)) for the chain geometry (6.3) 

TT 



Here we use the results of Bariev and McCoy et. al. [B2l p3J, and along the 
defect have changed the bond coupling from Kc ~ log(i(l -I- v^)) to Kc + A. 
This is not quite correct, because A is the coefficient of the energy operator in 
the continuum formulation, while K = Kc + A is the coupling in the lattice 
formulation. However, taking this difference into account will only give a (A- 
dependent) rescaling of our renormalizaion group flows, which will not affect the 
qualitative results. Some subtleties regarding the branch of the arctangent in 
Eq. (^.2|) for antiferromagnetic ladder couplings are dealt with in Appendix O. 



We get the OPE coefficient fecrcre = —5(A) from |34 



a{x - ^-)a{x + ^) = 5-s' [l + 5ge+.. .] (6.4) 

We can now get the renormalization group equations to 1-loop order solely from 
the OPE's pll, Isl : 
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^ = (l-.g(Af)A-8A2 (6.5) 

^ = -4.9(A)A^ (6.6) 

The flows for the ladder and chain cases are shown in Figures 5 and 6. The 
flows for the ladder case show that perturbations about the point A = A = 
(the defect-free point) eventually flow to the point with (A, A) = {0,—Kc)- 
This value of A corresponds to vanishing bond strength along the ladder. So the 
renornialization group flow takes us to a point with two decoupled half-plane 
Ising models and no random magnetic field. 

We can check our result by looking at the flows around the decoupled 
point. At this point the spins can be represented by free fermions, and we 
can calculate the renormalization group equations about this point. These re- 



sults agree with the flows around the decoupled point obtained from Eq.(6.5) 
and Eq.( |6.6[ ). In making the comparison, it is important to note that in the 
replica formalism, disorder in the magnetic field corresponds to an operator 
Sq=i X]/3=i '^a{x, 0)(T;3(a;, 0), but in Eq.(6.1), A is the coefficient of this opera- 



tor with the a = P terms removed. This means that A is not really the strength 
of the random magnetic field, but a linear combination of the strength of the 
random magnetic field and the bond stength (A). So, in figure 5, adding a ran- 
dom magnetic field to the decoupled point, (A, A) = {0,—Kc), moves us to a 
point with A > and A > —Kc, and the renormalization group flows from this 
new point eventually go back to the original decoupled point. 

The flows in the chain case take perturbations about the defect-free point 
to (A, A) — (0,— oo). Again, the random magnetic field has vanished at the 
new fixed point. To interpret this value of A, we look at the spin opera- 
tor. The dimension of the spin operator at A = — oo is given by Eq.( |6.3[ ) as 
Xcr — ^(?(— oo)^ = i, which is the same as the dimension of the spin operator 
along the edge of an Ising model with free boundary conditions. The possible 
boundary conditions of the Ising model have been completely classified [Q and 
the only boundary condition where the spin operator has dimension 1/2 is the 
free boundary condition. This makes sense physically, because A = — oo cor- 
responds to an infinitely antiferromagnetic coupling that produces alternating 
spins along the defect, which upon coarse-graining gives net magnetization zero 
everywhere along the line. We conclude that in the chain case we also flow to 
two decoupled Ising models with free boundary conditions. 

If is not hard to see that adding higher cumulant terms of the magnetic field 
to Eq. ( |6.l| ) will not change the qualitative results of our 1-loop renormalization 
group calculations. 

So far, we have represented the perturbation in the energy operator as chang- 
ing the bond strength either in the vertical or the horizontal direction, but 
not in both. More generally, we should represent the perturbation as chang- 
ing bond strengths in both directions. However, analagously to the ladder 
and chain cases, we expect that a more isotropic treatment would only give 
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a different monotonically decreasing function g{X), and that the point with 
(A, A) = {0,g~^{l)) would still be a stable fixed point with a large basin of 
attraction (including the defect-free model at A = A = 0). And by the classifi- 
cation of Ising model boundary states in [^, the point with (A, A) = (0, g~^{l)) 
will always correspond to two decoupled half-plane Ising models with free bound- 
ary conditions and no random magnetic field. 

7 Conclusions 

We have found a new universality class of disordered defect lines. The defect 
lines exist in various two-dimensional Statistical Mechanical models, such as the 
Tricritical Ising model and Tricritical 3-state Potts model. The large- distance 
behavior of these defect lines has been shown to be disorder-dominated. Two- 
point functions along the defect exhibit multifractal behavior, and universal 
(normalized) amplitudes of one-point functions are non-self-averaging. We have 
also argued that when a random magnetic field is applied along a single line of 
the Ising model, it causes the two sides of the defect to decouple, and to turn 
into two half-plane Ising models with free boundary conditions and no disorder. 

Results for the defect line in the m}^ Virasoro minimal model were obtained 
by a 1/to expansion. However, the physically most interesting models are at 
low m. For example, m = 4 corresponds to the Tricritical Ising model, with 
a random bond strength (or a random chemical potential) along a line. Our 
calculations show that this results in disorder-dominated long-distance behavior. 
It would be interesting to understand this model in a more fundamental and non- 
perturbative fashion. The random boundary/defect fixed points that we have 
found in this paper are, besides the bulk random bond g-state Potts models, a 
rare case where detailed analytic information about random critical behavior is 
available. In particular, it would be most interesting to compare our results for 
the random defect lines in minimal models with future numerical results, such 
as those obtained by Jacobsen and Cardy for the bulk random Potts modelsp2|. 

As a byproduct of our calculations for the defect line, we have also obtained 
multifractal energy-energy correlations in the bulk random bond q-state Potts 
model. 

Finally, we comment on the R.G. analysis performed under the assumption 
of broken replica symmetry. Such a calculation was done in |3^ for the q- 
state Potts model with bulk disorder in the bond strength, where it was found 
that the replica symmetric disordered fixed point is unstable to a new fixed 
point with broken replica symmetry. However, numerical tests show that the 
Potts model with random bonds is best described by the replica symmetric 
fixed point |g^, ^. An identical analysis to that of g^ for the random defect 
problem that is the subject of the present paper shows that, again, the fixed 
point considered in this paper is unstable, and flows to a new stable fixed point 
with broken replica symmetry. Thus, our random defect problem may provide 
further insights into the significance of the replica broken fixed point. 
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A Renormalization of A 

We calculate the renormalization of A to second order. The structure of the 
calculations closely parallels that in [p^. We get the coefficients in 

A(r) = r^ (Ao + A2(r, e)Ag + (A3i(r, e) + A32(r, e) + ^33(r, e))Ag + . . .) 

(A.l) 
where, as explained below, the different terms come from the different types of 
contractions at each order of perturbation theory. All integrals are regularized 
at short distances by analytic continuation in e = — q-r and at large distances 
by an infrared cutoff r. 

A.l First Order 

To lowest order we bring two perturbation terms together 

.2 



A2 !■ 

-f J dx,dx2Y.<l)U^l)<l)l2i^l)J2f,,ix2)<l)U^2) (A.2) 

We get the same perturbation back again (and thus a co ntribution to A) when 



/3 = 7 and a ^ S. This gives us the A2 term in Eq.( [A.f| ) 



A2 = 2{n^2) / dx2 < <l)i2{xi)(l>i2{x2) >o= Hn - 2)^ (A.3) 

\3:i—X2\<r 

The subscript on the correlator indicates that it is calculated in the defect-free 
theory. 

A.2 Second Order 

To second order we have 



2_^<Pi2[X2)(pi2(X2) 
7#5 



— - / / dxidx2dxs 

|a:i— a:2t<^ \xi—X3\<r 


n 
)_^C^U^l)cf>f2i^l) 







(A.4) 



This gives us several possible contractions. We get one possible contraction 
when (3 — "/, S = 11, a ^ S, f3 ^ i>, a j^ v: 



A3i=4(n-2)(n-3)Ji, (A.5) 
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where we have defined 

Jl = / dx2 / dxs < (t)i2{xi)(f)i2{x2) >0< 4'12{X2)4'12{X3) >0 

\x2—xi\<r 12:3— 2:i|<r 

= dx2 dX3 \ Xi - X2 \~^'^^\ X2 - X^ \~'^~^'' 

\3:2-'Xi\<r |£i;3— £Ci|<r 

= 2r^' f dy\y\'^+^'' f dz \ z \-'^+'\ I - z \-^+'' (A.6) 

Jo "'|z|<i/y 

We have transformed to coordinates y — ^^~^'- and z ~ £3;:£i \\/q can now 

^ r X2—X1 

extend the integral over z to go from —00 to 00, since this will only change Xi 
to 0(6°). This gives us two integrals which we can do exactly : 

= 4^+O(e0) (A.7) 

Note that in the integrals we have picked out one point xi as special, and 
integrated over X2 and X3, but in the combinatorial factor 4(n — 2)(n — 3) we 
have treated all points symmetrically (i.e., have treated a — S, P — fi the 
same as /? = 7, (5 = /i). This is permissible to this order in perturbation 
theory, since we are only concerned with the poles in e, which result from short 
range divergences and are the same for any permutation of the contractions. 
We can see this explicitly by considering a different arrangements of the same 
contractions: 



X 1 = I dx2 I dxz < (J)i2{xi)(t)i2{x2) >0< (t)12{xi)(j)l2{x3,) >Q 

\x2—xi\<r |a;3— a;i|<r 

= 4^=Xi+O(e0) (A.8) 

We get a second possible set of contractions when a — j = fi and 5 — v ^ j3: 

A32 = 4(n - 2)X2 , (A.9) 

where we have defined 



I2 = / dx2 i dxs.[< (t>l2{xi)(j)l2{x2)(l)l2{xs,)(l)l2{Qo) >Q X 

\x2-xi_\<r |a;3-£Ci|<r 

< (I)l2{x2)cj)l2{xz) >0 -(< (J)12{X2)'J)12{XZ) >of] (A.IO) 
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The four-point function gives the coefficient of three (/)i2's projecting to a single 
012. We have subtracted off (< ^12 (2^2)012(3^3) >o)^ - this term corresponds 
to the contribution from two perturbation terms (disorder operators) getting 
close, and does not affect the dimension of the operator at xi. This term is 
only a contribution to the free energy, so subtracting it off simply corresponds 
to normalizing the correlation functions. We get the four-point function from 
the Coulomb gas formalism Eol Ell: 



< 012(O)012(l)012(z)012(^) >0 = | Z r'+1 1 - Z T"'/' 



^ I z |2-4'/3| F(l - 1, 2 - 6; 2 - -; z) |2 -f I i^(l - 1, 1; -; z) |2 
4 ' ' ' ^ 3' 3 ' '^' ' ^ 3'3' 3 ' ^' 



,(A.ll) 



where we have defined 



fci = 



4(r(f))^r(2~6)r(i-f) 
3r(i)r(-i + 6)r(2-f) 



= l + 0(e) 



(A.12) 



and F is the hypergeometric function. If we take the limit as e — > at fixed z 
for the four-point function, we get 



< 012(O)012(l)012(2;)012(oo) >0^^ 

sign(z) sign(z — 1) 



sign(z)sign(z — 1) 



|z-l 



z\ 



(A.13) 



We now go back to X2, which we can calculate by using the symmetries 
X2 — > —X2 and X2 <-> 2^3 to cut the integration region down by a fourth, and 
then transforming to new coordinates y — ^^~^^ and z = ^^~^^ . The integral 
now exactly factorizes into two one-dimensional integrals : 



X2 = 4r 



1 /.I 

2£ / j„, „,-l+2e 



dy y~ 



dz [< 012(0)012(1)012(2)012(00) >o X 



< 012(1)012(2) >0 -(< 012(1)012(2) >o)'] (A.14) 



The y-integral gives ^ , while the z- integral can be rewritten as 



dz 



< 012(0)012(1)012(2)012(00) >o< 012(1)012(2) >o 



- < 012(0)012(2) >o - (< 012(1)012(2) >o) 



+ I dz < 012(0)012(2) >o 



(A.15) 
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The second integral is exactly -. It is straightforward to use Eq.(A.ll) to 
check that as z ^ 1, the four-point function goes to | z — 1 |^i+' +0{z — 1). 
So the first integral converges everywhere, and we can get it's value to 0{e'^) by 
simply replacing e with zero inside the integral: 



dz 



1 / . / X . / -,x sign(z) sign(z— 1), 
(sign(z)sign(z - 1) + , \ ^| - ^ '-) 



z~\ 



z-\ 



^1 



\z\ |z-ip 
-l + 0(e) 



0(e) 



(A.16) 



Putting this all together gives 



J, = 4r^ 



-2r' 



2e 



1 + 0{e) 



1 2 



(A.17) 



Again, as with A^i, other possible permutations of this contraction like 



J2 = / dx2 / dx^i < (i)i2{xi)(l)i2{x2)4>12{x-i)(t)i2{Qo) >o 

\x2—xi\<r |a:3 — a;i|<r 

< (I)l2{xi)(t)l2{x2) >0 (A.18) 

give the same result up to terms of 0{e^). 

Finally, the last possible contraction in Eq.(A.4) comes from a = 7 = /i and 
(3 = 5 = v. This is 



A 



33 



dx2 I dx3 

\x2—xi\<r |x3 — a;i|<r 



(< 0i2(a;i)(/)i2(a;2)0i2(a;3)0i2(oo) >o)^ 



(< </'i2(a;2)0i2(a;3) >o)' 



16 2e /■' dy /■! 



y J-1 



dz 



(< '/'12 (0)012 (1)012 (^;)<?!>12 (00) >o) 



(< 012(1)012(2:) >o)^ 



(A.19) 



The y-integral gives |. The z-integral is evaluated similarly to Eq.(A.15). We 
rewrite it as 



dz 



(< 012(0)012(1)012(^)012(00) >o)' 
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dz 



- (< ^12i0)<pl2{z) >of - (< (/'12(1)'/'12(2) >of - 1 



l + (<'/'l2(0)<^12(^)>0)' 



(A.20) 



The second integral can be evaluated exactly, and is 0{e). The first integral 
is nowhere divergent, so we can get its value to 0{e'^) by taking e to inside 
the integral. Then squaring Eq.( A.l^ ) gives 



(< 012 (0)(/)12 (1)012 (2)(/'12(oo) >o) 

e^o f . . ^ . . -,^ , sign(z) sign(z-l) 

sign z)sign z - 1) + — . . 

z- 1 \ z \ 



= 1 + ^ 



(z-1)^ 



(A.21) 



So both z-integrals are 0{e), and ^33 is C(e°)- Since A33 has no pole in e, 
it can be dropped to this order. Combining all these results gives Eq.(3.2). 



20 



B Renormalization of A 



NMn 



To get the dimensions of the moments of (/)i2, we need the coefficients of Aq 
and Al in ZNMnir), where ANMnir) = ZNMn{r)ANMnfi- The integrals in 
this section wiU be the same as the integrals in appendix |A[ The only difference 
will be in the combinatorial factors which precede them. In this section we will 
use {N, M, n\ and \N, M,n) to represent the bra and ket forms of the operator 



Onmu defined in Eq.(4.2). However, before we start, we should make sure 
that no problems occur with mixing, either for the defect line considered in 
this paper, or for the parallel calculation in the q-state Potts model with bulk 
disorder in the bond strength. 

B.l Mixing With Other Operators 

The combinations Onmu of Eq.([4.2|) diagonalize the operators of the form 
]^(/>"2- However, we also need to check that these operators don't mix with 
other operators that have the same scaling dimensions at 771 = 00. We only 
need to consider operators of the form 0ig, because these operators form a 
closed subalgebra in the minimal models. To 2 loops, only 4 (^12 operators can 
be affected by the perturbation, and 4 dim((/)i2) = 2 at m = 00. It easy to see 
that the only operator, or descendent of an operator, or combination of opera- 
tors, in the 4>iq subalgebra with dimension 2, is (j>i3 (therefore, to this order, no 
mixing occurs with derivative operators). (^13 does in fact mix with (0i2)^, and 
so (f>i3{4>i2)^~'^ mixes with (0i2)^. However, the 2-loop overlap integral is 



dX2 / dX3 < 4>12{xi)4>i2{x2)4>13ioo) >0 

|a;2— a;i|<r |a:3— £ci|<r 

< (f>12{xi)(l)i2{x2) >0 (< (t>12{xi)(l)i2{x3) >af (B.l) 

This has no poles in e, so we would only need to take this mixing into account 
to do 3-loop calculations. 

Mixing is potentially more of a problem when we use these calculations to 
get moments of the energy operator in the q-state Potts model with bulk bond 
disorder - this model is analyzed by expanding about the Ising Model |g^, ^ , 
so we need to look for mixing with other operators that have the same scaling 
dimension at q = 2. Denoting the energy operator in the a replica by e", we see 
that e°'e^e'^ and V^e" have the same dimension at q = 2, and that in general 
we can replace any two e's with a V^ without changing the dimension. These 
operators mix to two loops, with overlap integrals such as 



21 



(Px2cfx3 < e(0)e(x2)e(x3)V^e(oo) >o 

< e{0)e{x2) >o< e{0)e{x3) >o (B.2) 



This means at the combinations Onmu of Eg. ([4.21) (with 0i2 replaced by e) 
are no longer generally multiplicatively renormalizable, but mix with derivatives 
oiO j^_2 ]yj „• Luckily, however, no problem with mixing occurs for the operators 
ON,N,n and ON^N-i,n, which we know from the one-loop calculations to provide 
the leading and subleading decay exponents. No problem with mixing occurs 
for On N,n because the overlap 



^ e"e'' ^ e^e^\N, N, n >^ V^\N - 2, Af , n > (B.3) 

is for all M, where the V^ is understood as acting on only one of the e's. It 
is not hard to see this by counting the connected contractions of 



<N,N,n\Y^ e°'e^ ^ e^e* |iV -2,M,n> (B.4) 

If we expand | A'^ — 2, M, n > into monomials, each monomial term gives a con- 
traction of ~ this is because after contracting all the e's in the monomial 
term, we will be left with at least two antisymmetric terms from the < N, N, n\, 
and thus positive and negative contractions will cancel. So the contractions 
give a total of (even before the replica limit n ^ is taken). This is true 
even if the |A^ — 2, Af, n > has a V'^ on it, and even if < A^, A^, n| is replaced 
with < TV, TV — 1, n|. So mixing is not a problem for the ON,N,n and Ojq.N-i.n 
operators in the q-state Potts model with bulk bond disorder. 

B.2 First Order 

Returning to our defect model, the 1^* order correction comes from 

Y. ^i^^i^ I 1^' ^^' ") -^ ~^^Mn \N, M, n) (B.5) 

where we contract one 0i2 operator in the Y^^j^a <^?2<^i2 with one (j)i2 operator 



in the |iV, Af, n). The integral is the same integral done in Eq.(A.3), giving 2—. 
We can get the combinatorial factor bNMn by contracting with (A'^, Af, ri| on 
both sides : 
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^^*'^" - {N,M,n\M,N,n) ^^-^^ 



This combinatorial factor is found in subappcndix IC.4 Leaving our result in 
terms of bjqMn , the order Aq term in Z^mu is 

2bNMn (/-^\ Ao (B.7) 

B.3 Second Order 

To second order, we need to look at the number of ways that we can take 

n n 

< N, M, n\,J2 '^?2'^?2 E '^?2'/'i2 ^< N, M, n |, (B.8) 

The index /i in < iV, M, n j^ is a label used to make discussion easier, and does 
not signify an independent variable. The ^ simply indicates that for discussing 
contractions, the dummy indices Ui in Eq.(4.2) will all be called /ij for i = 
l,2...n. 

There are now a number of possible contractions to consider. The one with 
Hi = a = ^ = Vj and ^k = P = S = ui (for some i,j,k,l) gives an integral over 
a squared four-point function. We did this integral for A33 in appendix |^ and 
found that it had no singularities. So we can ignore this term. 

We now consider the contraction with /i^ = a = 7 = za, and P — S, which 

(2) 
gives a combinatorial factor that we call Qf^Mm ^^"^ ^^'^ contraction with jii = 

P = "f = Vj, Hk = ct, and 5 = vi, which gives a combinatorial factor that we call 
Qnmu- ^^ ^^ ^'^ ^^ understood that no replica indices other than the specified 
ones are equal to each other. In both cases, we get the same integral as in 



Eq.(|Al0|) (or Eq.(|Al8|)), which was found to be I2 = -2r2^(i - 4)- 



To get the combinatorial factor for Q'j^j^j^^ we look at 

^< JV, M, n I Y.l^p 0?,0f, YT,^, 4>l2^i^ I N, M, n >. 

< TV, A-f , n I A^, M, n > ^ ' ' 

and only allow the contractions described above. For the contractions allowed 

f2) 

in Qj^Mn^ ^^^ replica indices that appear in a monomial term must be the same 
in the bra and ket sides. The a = 7 replica index can be the same as any of the 
these A'' monomial terms. The P — 5 replica index must be something different, 
so has [n — N) choices. We also get an overall factor of 4 for the different pairs 

of terms that we could have chosen to contract. So Q)^j^.jj^ = 4iV(ri — N). 

(3) 
Qnmu *^^^ ^^ calculated similarly. The a index contracts to the left and the 

S index contracts to the right - this gives a combinatorial factor < A^, M, n \ 

T,2^5 </'i2</'i2 \ N,M,n>= Inmu < N,M,n\ N,M,n>The P = -f index can 

then be any of the (TV — 1) replica indices in the right (or left) | N, M, n > 
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(excluding the two which are equal to a or f3). And there is again an overall 

(3) ~ 

factor of 4, so Q]vm« = 4(iV - l)bNMn- 

Now consider the contraction with fii = a, (3 = j and S — Vj, which gives a 

combinatorial factor Qj^jij^i ^-^d the contraction with ^ — a, /ij — j, i/i = f3, 
and i^j = S, which gives a combinatorial factor Q^j^.j^- Both of these terms give 
the same integral as in Eq.( |A.6| ) (or Eq.( [A.8|) ), and thus a factor of Xi = 4^ . 
Qnaiu '^^^ ^^ evaluated in the same manner as Qj^^^^ S^'^^^S Qnaiu ~ 
4(n — N — l)b]\[Mn- The term Q^vMn i^ more complicated and is found in 

subsection p. 3| to be QpfMn — ~4A^(n — iV) — 2(n — 2)feAr7Vfri+ (^AfMn)^- Putting 
this all together, we get the result 



ZnMu = 1 + Z^ATAfn — Ao + -^{QpfMn + Q N Mn)i^'^ ^ ^^") 



^(o^L+Q^L)(4j)+o(A^) 



= 1 + 26wM„^Ao - 4 (7V(n - iV) + (A^ - l)6ArM„) ^Ag 

+ (2(6jvM„)' + 4(n - 2)bNMn) ^A2 + 0(A3) (B.IO) 
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C Combinatorial Factors 

In this appendix we calculate the combinatorial factors used in appendix H. We 
want the number of possible contractions between the irreducible representations 



of Sn, given by Onmu in Eq.([4.2|), and copies of the disorder operator. The 



trivial spatial dependence is supressed in the equations below. 

C.l Normalization 

First, we calculate the normalization of | iV, M, n > 

A[N, M, n] = <N, M, n \ N, M, n >= ^ ^ 

l<Qi<(n-M) l<ft<(n-M) 
(0?^ - r^) • ■ • (0?!' - 012"^''"'V?I'"^ • • ■ <l>'l2 > (C.l) 

If we expand out the two terms {(j)"" — '/>i2~ ~ ) and {4>{2 — <th2 ) i^ito 

monomials, the cross terms are 0, and the remaining terms are of the same form 
A[. . .] as before, but with different values of M, N and n. We get 



A[N, M, n] = A[N, M-l,n-l] + {n~M~N+ lfA[N-l, M- 1, n- 1] (C.2) 
We can calculate A[N, M, n] explicitly when M = : 

A[N,0,n] 



E E <<2<2- 


•.'^?2" 


\ 012 012 • 


■42> 


l<ai<nl<ft<n 









N 

= 7V!j|(n-iV + i) (C.3) 

i=l 

Here the N\ comes from the number of ways to contract the left side with the 
right side, and the factors of n — A^ + z come from the different ways to pick the 
Ui once the contractions have been done. Given the value of A[N^M,n] when 
M = 0, and the recursion relation above, we can show by induction that the 
general solution is 

N M 

A[N,M,n] = {N - M)\\W{n - M - N + i)]\W{n - 2M + I + i)] (C.4) 
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C.2 OPE - 1st order term 

Here we want to calculate 



B[N, M, n] EE< TV, M, n I y 05*2012 \N,M,n> (C.5) 



/ , V^12V^12 



Note that three-point functions of (/)i2 are by the (/)i2 -^ — 0i2 duality 

6?f - r,2^^'-'^) and {4^ - 



symmetry. As before, we expand out the {4>l2 ~ 4>i2 ) and '^^'^'^^^ 



•5^12 ) terms into monomials. The direct terms again give contractions 

of the form i?[. . .], but we now have a cross term in which either a or /3 must 
be equal to n — (Af — 1). We get 



\b[n,mm = 


iB[iV,M-l,n- 


-1] 


+ (n-M- 


N + lf-B[N - 


1,M- l,n- 1] 


-2(n - M 


-N + 1) >; 


n 

>: ^ 




l<a,<(n- 


-M) l</3.<(n-M) 


< {4>i^ 


-ri2)-{<t>i2'-'- 


-</.;'2"^*'''V?2"'/'r2"" 


K'/'S 


^^12)-{^fr'' 


-0;'2-(^-)) <^fr 



*?2" I '/'?2 



6^2" > (C.6) 

The contraction in the last term looks just like the contraction A[N, M — l, n— 1], 
with a taking the place of /3m- The only difference between the sum above 
and A[N,M — l,n — 1], is that the sum over a above includes terms with 
(n — M + 2) < a < n, and terms where a is the same as some other term Pi, 
for I < i < (M — 1). However, these two extra contributions are equal and 
opposite, being equal to ±(M - l){n - M - N + 1)A[N - l,M-2,n-2]. So 
they cancel, and we have 



-B[N,M,n] = -B[N,M -l,n-l] 



{n- 


- M - 


- N + lf-B[N - 


l,M - 


-l,n 


2{n 


-M 


-N +l)A[N,M 


- l,n- 


-1] 



(C.7) 



As with A[N, 0, n], it is easy to calculate B[N, 0, n] : 



N 



B[N, 0, n] = 2N{Nl)[[[{n - N + i)] = 2N{n - N)A[N, 0, n] (C.8) 

4=0 

Given the recursion relation for B[N,M,n], the initial condition B[N,0,n], and 
the result for A[N,M,n] in the previous section, we can show by induction that 
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C.3 2nd order term 

We want to calculate 

n n 

C[M, iV, n] = 4 < iV, M, n I ^ 0^20^2 ^ ^\^^\^ \ N, M, n > (CIO) 

where we require that a and 7 contract to the left, the /3 and (5 contract to the 
right, and a, /?, 7, (5 are all distinct from one another (other possible directions 
of contractions give the factor of 4 in front). This combinatorial factor arises 
in the contraction for Qf^j^j^ in subappendix |B.3| As in the previous subsec- 
tion, expanding out the (^"2" ~ '^12 ) ^^'^ {4'i2 ~ '^iT ) terms into 
monomials gives back two terms of the form C[. . .], and a more complicated 
cross-term. The cross term almost has the same form as B[N, M — 1, n — 1], but 
we also get some extra terms because the sums over a, /?, 7 and S are different 
than we would have in a -B[. . .] term. Counting all the ways in which our cross 
term differs from B[N^ Af — 1, n — 1] is tedious but straightforward, and we get 

-B[N,M -l,n-l] - {M - \){n~M - N + \fA[N -1,M ~2,n~2] 
-(M- l)yl[iV,Af-2,n-2]- {n - M - N)A[N,M -l,n ~1] 
^ -B[N, M -l,n-l]-{n-N - 1)A[N, M - 1, n - 1] (C.ll) 

The recursion relation is 

-C[N, M, n] = -C[N, M - 1, n - 1] + 

(n - M - TV + lf-C[N - 1, M - 1, n - 1] - A{n - M - N + 1) 

-B[N, M -l,n-l]-{n- N- l)A[N, M - 1, n - 1] I (C.12) 

Combined with the initial condition, 

C[7V, 0, n] = 4:N{N - l){n - N){n - N - 1)A[N, 0, n], (C.13) 

we can find the value of C[N, M, n] for all M and N by induction. The result is 

where bNMn was defined in Eq.(|C.9|). 
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D Ising Ladder Defect for K K 



The branch of the arctangent used in Eq.(3.2) for antiferromagnetic {K ~ Kc + 
A < 0) ladder couphngs requires some explanation. We take the value of the 
arctangent to be in (0, ^) if its argument is positive (i.e. K > 0) and to be in 
(^,7r) if its argument is negative (i.e. K < 0). This makes the slope of g{X) 
continuous through K = 0. On the other hand, the results of [^2[ |3^ have a 
slope discontinuity at K = 0, and have g symmetric under K — > —K. This 
slope discontinuity results from a level crossing in the lowest scaling dimension 
for operators on the boundary |Q . If we let at be the spin on one side of the 
defect, and ab the spin on the other side, we see that the operator with the 
lowest scaling dimension changes from at + at to at ~ ab as K goes through 0, 
so that while the dimension of each operator changes smoothly through K = 0, 
the dimension of the lowest scaling operator does not. However, if we take our 
random applied magnetic field to not vary across the defect, it couples to at+ab 
only, and we want the lowest scaling dimension for K > 0, but the 2"*^ lowest 
scaling dimension for K < 0. We can get these from Q, thus justifying the 
branches of arctangent chosen above. 

Note that if we had used the other branch of the arctangent, corresponding 
to a magnetic field uncorrelated across the defect, the flow picture in the ladder 
case would have been symmetric under K — > —K, and the entire (A, A) plane 
would have flowed into the decoupled point. 
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Figure Captions 

Fig. 1. A two-point correlation function for operators lying along the defect. 

Fig. 2. A one-point function for an operator in the bulk. 

Fig. 3. The ladder defect in the Ising model. 

Fig. 4. The chain defect in the Ising model. 

Fig. 5. Renormalization Group flows for the ladder defect in the Ising model. 
A is the strength of the disordered magnetic field along the defect line, and A 
is the bond strength along the line. 

Fig. 6. Renormalization Group flows for the chain defect in the Ising model. 
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